Excitation of Neutron Star Oscillations by an Orbiting Particle 
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The excitation of neutron stars is expected to be an important source of gravitational radiation. 
Of fundamental importance is then to investigate mechanisms that trigger oscillations in neutron 
stars in order to characterize the emitted radiation. We present results from a numerical study of 
the even-parity gravitational radiation generated from a particle orbiting a neutron star. We focus 
our investigation on those conditions on the orbital parameters that favor the excitation of ui-modes. 
We find that, for astrophysically realistic conditions, there is practically no w-mode contribution to 
the emitted radiation. Only for particles with ultra-relativistic orbital speeds > 0.9c, the w-mode 
does significantly contribute to the total emitted gravitational energy. We also stress the importance 
of setting consistent initial configurations. 



I. INTRODUCTION 



Neutron stars as a source of gravitational radiation have been under investigation for over 30 years, with the 
foundation developed by the work of Thorne and Campolattaro |Q] . The radiation emitted by an excited non-rotating 
neutron star basically consists of a superposition of the characteristic oscillatory modes of the star [|^^ . Those modes 
fall into two categories. One class are spacetime modes, such as the w-modes and the trapped modes for very compact 
star. The other class consist of fluid modes, such as the /-mode and the p-modes which also exist in Newtonian 
stars. Due to the emission of gravitational waves, both types of modes are damped, with the spacetime modes much 
, stronger than the fluid modes. 

Recent studies of the time evolution of these oscillatory modes [^D have shown that generic, but somewhat un- 
physical, initial data can excite the first w-mode, the /-mode, and a few p-modes. However, the strength of these 
O modes depends crucially on the particular choice of initial data. Conformally flat initial data, for instance, has the 
^ '■ effect of almost completely suppressing the excitation of w- modes ||^, whereas the /-mode is generally present. Work 
5-H . by Andersson and Kokkotas Q shows that, by extracting the frequencies and damping times of the first w-mode and 
the /-mode, one could in principle determine physical parameters of the neutron star such as mass, radius and even 
^ , the equation of state. Of course, the success in the determination of such physical parameters will depend on the 
particular modes present in the signal and the noise level. 
. , The issue of interest for the work in this paper is whether astrophysically motivated initial neutron star perturbations 
■ are able to excite the w-modes, so that accurate determination of the mode parameters could be possible. One of the 
■ ■ ' most likely events that may excite strong oscillations is the birth of the neutron star. However, a study of this process 
would require a rather complete understanding of the collapse of the progenitor star within the full nonlinear theory 
of general relativity. The alternative is to study astrophysical systems involving neutron star oscillations that can be 
treated using perturbation theory. A system that allows such an approach and, at the same time, could be viewed as 
an approximation to a realistic astrophysical system is that of a "small" particle orbiting a neutron star. Here small 
is to be understood in the sense that the mass /i and size of the particle are much smaller than the mass M and radius 
-R* of the neutron star. Under this approximation, the particle moves along a geodesic in the spacetime generated by 
the neutron star, causing metric perturbations that propagate and eventually excite the neutron star. Of course, in 
principle those gravitational perturbations in turn would cause deviations of the particle's trajectory; however, these 
back-reaction effects are of second order and will therefore be neglected. 

The specific goal of this paper is to investigate the gravitational radiation emitted from the orbiting and scattering 
of a particle by a neutron star. Collisions of particles with the neutron star ^ are not considered because it is not 
clear at all how to treat the impact and subsequent merge of the particle with the neutron star. Contrary to previous 
studies of excitations of neutron stars by orbiting particles, we perform our calculations in the time domain instead 
of the frequency domain jl^ In particular, we use the linearized ADM- formalism jF~ 
for the numerical initial value problem. The relevant equations have been derived in 
paper I. Considering the dynamics in the time domain has as a consequence that we are forced to "smooth out" the 
particle; that is, the 5-functions in the sources of the equations due to the presence of the particle are approximated 



, which is naturally adapted 
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by Gaussians. We show, however, that our model of the particle is self-consistent and convergent. Another important 
issue is the prescription of appropriate initial data that satisfy the constraints. 

The paper is organized as follows. In section ||, we present the equations to describe the unperturbed configuration 
and its perturbations, as well as the method to treat the orbiting particle as a perturbation of the background 
spacetime. Section III deals with the numerical approximation of the particle terms in the equations. In section |rv|, 
we outline the problems with choosing appropriate initial data, and, by referring to the flat space problem, we give an 
approximate prescription of "good" initial data. In sections ^ and VI, we present and discuss our numerical results. 
Details concerning the derivation the particle terms can be found in the appendix. 



II. FORMULATION OF THE PROBLEM 



The energy-momentum tensor T^'^ of a point particle with mass /i is given by 



T^"' =^l / ^f-*) (x'" - A''(r)) u^'u''dT 



where is the 4- velocity of the particle 



cose(t)) , 



dT ' 



(1) 



(2) 



and T the particle's proper time along its trajectory X'^{t). 

We now orient the coordinate system in such a way that the particle's orbit coincides with the equatorial plain 
of the neutron star (8 — 7r/2). Besides, the particle's coordinate time T is identical with the time coordinate 
t of the spacetime in which it moves. Therefore we will use t to parametrize the path of the particle A'^(i) = 
[t, i?(t), n/2, 

From the geodesic equations '^=^j— = U^U^.j^ = we find: 
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(3a) 
(3b) 
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where E and L are the energy and angular momentum per unit mass of the particle, respectively. We also recall that 
for the Schwarzschild metric we have 



e^-' = e-^^ = 1- 



2M 



We can use (^^ to eliminate the proper time r from equations (|^) and (|^): 
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Eqs. (|5^) and ( pb| ) can be used to replace the quantities Vr = ^ and t;^ = ^ in the source terms of the particle. 
However, we also have to explicitly solve them for the particle's trajectory coordinates R{t) and ^{t), since we need 
those coordinates in the J-functions 6{r — R{t)) and S{(p — ^{t))- 

The background spacetime is solely due to the neutron star, which we model as a non-rotating polytrope with F = 2 
and K — 100 km^. Outside the star, the spacetime is given by the Schwarzschild metric. The equations governing the 
stellar oscillations are derived in paper I. After expansion in Regge- Wheeler tensor harmonics, these equations can be 
grouped into two families, one describing the axial perturbations and the other describing the polar perturbations. 
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In the exterior region, we now have to incorporate the source terms due to the particle into the perturbation 
equations. For that purpose, the energy-momentum tensor of the particle has to be expanded in tensor harmonics. 
The derivation of this expansion is given in the Appendix. There is a slight inconvenience created by the fact that the 
perturbation equations are expanded using the Regge- Wheeler harmonics which do not form an orthonormal set, 
whereas the energy-momentum tensor of the particle is expanded using the orthonormal harmonics given by Zerilli 
Ip^ . However, the conversion is straightforward. 

Even though the main focus will be on the polar perturbations, we also present the axial equations since they are 
so simple. In the Regge- Wheeler gauge, the only non-vanishing axial perturbations of the shift vector and spatial 
metric are respectively given by 
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Similarly, the extrinsic curvature perturbations read 
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with Wim and Xim defined in Eqs. ( |A6| ) and (A5) of the Appendix. With the inclusion of the particle terms the 
evolution equations read (see Paper I) 
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The momentum constraint relates the extrinsic curvature coefficients to the particle's source term via 
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For completeness, we also list the source terms of the particle: 
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tj — e 
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(12) 
(13) 
(14) 



For radial infall, L = 0, and all source terms vanish. Hence, in this case the radiation is of pure even parity. 

Let us now turn to the polar perturbations. Here, we have to use some caution in including the particle terms. From 
the ADM split, we see that the matter terms enter only in the equations for the extrinsic curvature coefficients Ki. In 
the Regge- Wheeler gauge, we are only left with equations for the coefficients Ki, K2 and if 5. In paper I, we showed 
that we have to choose initial data with K^^ — 0. Without the particle, the evolution equations then guaranteed the 
vanishing of for all times. However, in the presence of the particle, this is not true anymore since satisfies the 
following evolution equation: 



dt 



(15) 
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This means that during the evolution K4 will be become nonzero even if it was set to zero initially. A non-vanishing 
K4 implies nonzero metric components Ti and V3; this violates the Regge- Wheeler gauge condition. A solution is to 
choose a different lapse function perturbation a. If we pick 



1 



a 



T 



rS + leTTtg , 



(16) 



the last term exactly cancels the right-hand side of ( p^ ) and the vanishing of A'4 and therewith Ti and V3 is, again, 
guaranteed. 

The remaining metric quantities are expanded in the same way as in paper I: 



= {e^^K2, 0, 0) Yir, 
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Or sin^ T 

and the expansion of the extrinsic curvature reads 
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In paper I, we have shown that it is possible to write the evolution equations as a set of two coupled wave equations 
only in terms of the metric components S and T. With the inclusion of the particle terms those evolution equations 
read in the exterior region 
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where the Pi and P2 denote the source terms of the particle 



Pi^t''-2t' + {5iy'--]t'+t5-2{3v' 
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with the U given in the Appendix. The evolution equations which govern the oscillations in the interior of the neutron 
star are given in paper I (e.g. (49) and (58)). Furthermore, the perturbations in the exterior are completely described 
by the Zerilli equation (eqs. (61) of paper I): 



Y{r)Z + Sz , 



with the following explicit form of the source term 

4iy M 
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where 2n = l{l + 1) — 2. In addition, we have three constraint equations, namely the Hamiltonian constraint, which 
now reads 



T" + v'T' -rS'-[b- + e2Ai([±l) 



1 
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T-[2rv' + 2 + -e'H{l + l)]S ^ -Sne^^^— S{r ^ R{t))Yi*^ , (26) 
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and two momentum constraints 



rK'^ - ie^^ 1(1 + l)K2- r^K - {rv' + l)K^ = Sne^^'^iEvr S{r - R{t))Yila (27) 



rK'^ -r'K-2K,^ 16^e2''-i!^^(, _ i?(t))y^;^ . (28) 



Notice that the quantity K2 appears in Eqs. (|2^) and (|2^) but not in the evolution equations. Therefore, we will 
eliminate it by differentiating ( [27| ) with respect to r and using ( p7| ) and (^8|) to eliminate K2 and i^Tg. The resulting 
equation is then second order in and reads 

^^+e-^W.-f2..' + 2+i.-' 
r / V 2 



ir^' + ly'K'^ ~rK'~[5-+ e^^ ' ' ] - [ 2rv' + 2 + -e^H{l + 1)]K 

(29) 



= Stt^ {e^^VrE5'{r - i?(t)) + imL5{r - i?(t))) y,;;, . 

It should be noted, that this equation also follows from taking the time derivative of the Hamiltonian constraint (|2^ ) 
and taking into account that S — K and T — K^. 

As usual, the initial data must satisfy the constraint equations. Unfortunately, there is no unique way to solve 
those equations. This is due to the fact that to a particular solution of the inhomogeneous equations, we can always 
add a solution of the homogeneous equation, which would correspond to adding some arbitrary gravitational waves. 
The problem of finding the "right" initial data that represent only the perturbations which are due to the presence 



of the particle and which do not contain any additional radiation will be discussed in more detail in section IV . 

For radial infall of aparticle from rest, Vr{t = Q) = Q and L = 0, and therefore t2 ~ H ^ 0. Thus the momentum 
constraints ( p7| ) and ( p^ ) can be trivially satisfied by setting the extrinsic curvature variables to zero. This situation 
corresponds to time symmetric initial data. We then are left with solving the Hamiltonian constraint ([2^). Of course, 
a particle falling from rest would fall radially towards the neutron star and eventually hit its surface. Since we want 
to avoid such an impact, we have to give the particle some angular momentum. In addition, we want to consider 
different initial radial velocities which means that we have to solve also the momentum constraints (|27|) and 



III. NUMERICAL IMPLEMENTATION OF THE PARTICLE 



The presence of a particle in the calculations introduces source terms in the equations that require the explicit 
forms of the spherical harmonics Yhn- The perturbation equations without particle are degenerate with respect to m 
since the background metric is spherically symmetric. However, the particle breaks this symmetry, and one is forced 
to consider the various m-cases. Fortunately, we do not have to consider all possible values of m for a given values of 
I since for negative m the spherical harmonics just undergo sign change and phase shift (Yj*^ ~ (— l)™y;^_m). The 
advantage of putting the particle in the equatorial plane (0 = ^) is that, in the even parity case, we only have to deal 
with multipoles with m — I, I — 2, ... , the remaining ones (m — I — 1, I — S, ...) have odd parity. Since the evolution 
code handles only real valued perturbations, we have to treat the real and imaginary parts of the spherical harmonics 
separately. Finally, all the equations will be solved on a finite grid, hence we have to approximate the (5-function by 
a narrow Gaussian 

S{r - R{t)) « e~ ^''^2^^ a small . (30) 

(TV27r 

To demonstrate the validity of this approximation, we have to ensure the convergence of the solution for u — s- 0. 
This can be done in two different ways. One can look at the convergence of the waveforms that are obtained in the 
evolution, or also one can monitor the violation of the constraints. A possible way do to so is to monitor the following 
quantity 

/ = / re^" (Ihs 

where the domain of integration is the region outside the neutron star. In the limit — s- 0, one gets that 1=1 
throughout the whole evolution. Numerically, we cannot take this limit with a fixed grid size since, eventually, we 
cannot sufficiently resolve the Gaussian. To check the convergence, we decrease both cr and Aa; by keeping the ratio 
a/ Ax constant throughout the sequence. By consecutively doubling the resolution and halving cr, we find that / — 1 



26)) dr 



(31) 
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approaches zero with second order convergence. For the numerical evolutions, we will use a value of (t/Ax = 0.15, 
which gives provides a suitable resolution of the Gaussian and its derivatives. 

There is still a subtle point. In deriving the source terms (see Appendix), we have tacitly transformed the particle 
coordinate R into the spacetime coordinate r since the presence of the (5-function makes a distinction unnecessary. 
However, in the evolution equations we have to take derivatives of the source terms with respect to r but not with 
respect to R, and therefore we would obtain different source terms if we had not changed the R's into r's. As an 
example consider the following two source terms 

Si{r)^ f{r)S{r-R) (32) 

and 

S2{r)^ f{R)6{r-R) , (33) 

which are equivalent because of the presence of the (5-function. But if we now differentiate Si and 5*2 with respect to 
r we obtain for 5*2 just the derivative of the (5-function, whereas for Si we also have to differentiate /. Analytically 
this does not make a difference, but if we approximate the (5-function by a Gaussian, then the two expressions for 
Si and S2 and their respective derivatives are different. To gain accuracy, we should have kept r'^dcp/drS^r — R) as 
Lr^ /R^6{r — R) and not just as LS{r — R). However, for the numerical evolutions the actual difference is negligible, 
so we have assumed the source terms to be of the form of Si. 



IV. SETTING UP THE INITIAL CONDITIONS 



Any construction of initial data for a particle initially located at (ro, $0) with initial radial velocity Vr and angular 
momentum L involves solving the Hamiltonian constraint (^). Usually, we set the initial angle $0 to zero. However, 
Eq. (^) contains two quantities S and T. This means that there is some freedom in choosing the initial values. We 
can e.g. either set the quantity T to zero and solve for S, or do it the other way round and set S — and solve for 
T. In the former case, we would have to solve a first order equation for S, and the latter case would lead to a second 
order equation for T. In Fig. |^ we show both possible sets of initial data. 

To assess which choice is the more natural, we consider a particle initially at rest in flat spacetime. Of course, the 
particle will remain at rest since there is no matter around that could attract the particle. Thus, the perturbation 
of the spacetime that is created by the particle will be stationary. Hence the equations of motion for the metric 
perturbations S and T will read ^ = and ^ = 0. From these conditions and from the Hamiltonian constraint, it 
then follows that S has to vanish. 

Of course, in the presence of the neutron star, those arguments do not hold any more, and S = will not be the 
right choice of initial data, but if the particle is initially far enough away from the neutron star, the error in setting 
S = should be very small. This error actually corresponds to an introduction of an extra amount of gravitational 
radiation that is not at all related to the radiation that is emitted when the particle moves through the spacetime. 
This extra amount will start to propagate during the evolution and eventually hit the neutron star and thereby cause 
it to oscillate. However, if we put the particle far enough away from the neutron star, the strength of the induced 
oscillation should be small compared to the ones excited when the particle comes close to the neutron star. 

In Figs. H and ^, we show the evolution of the two possible choices of initial data for a particle falling from rest. 
The left panels show the evolution of S and T in the case To = 0, whereas the right panels show the evolution of S 
and T with 5*0 = 0. The differences are obvious. The initial shape of T in the latter case is almost unchanged during 
the evolution, whereas S starts to acquire its "right" shape. In the other case, we can see a huge burst of radiation 
propagating in both directions. In the same time T is acquiring its "right" shape. In Fig. |^, we show the evolution of 
the Hamiltonian constraint where we evaluate the left-hand side of (|2^) which monitors the "path" of the particle. In 
both cases the graphs agree. Lastly, in Fig. |[ we also show the metric functions S and T at the end of the evolution. 
This figure clearly show that, regardless of the chosen initial data, the functions S and T will adjust to their proper 
values after having radiated away the superfluous initial wave content. 

We should note that we cannot escape the whole ambiguity of how to choose the variables S and T by using 
the Zerilli formalism, instead. There, any regular initial data are valid initial data since by construction the Zcrilli 
function always satisfies the constraint equations. Hence, we cannot see a priori whether or not the chosen initial 
data will have additional radiation content. 

If the particle is initially not at rest, in addition to solving the Hamiltonian constraint (|2^), we also have to solve 
the momentum constraints (|2^) and (|2^). Of course, as with the Hamiltonian constraint, we are again faced with the 
same kind of ambiguities in solving that equation. 
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Now, for a particle being initially at rest or very slow, or for circular orbits, we do not worry about what kind of 
initial data we choose. Since any initial gravitational wave pulse travels with the speed of light and is therefore much 
faster than the particle, it will long be gone when the particle comes close to the star. However, if the particle's initial 
velocity is close to the speed of light, then the particle "rides" on its own wave pulse, and it will be not clear any 
more whether the excitation of some particular modes of the neutron star is due the the particle itself or due to the 
initial burst which comes from inappropriate initial data. If the particle is slow enough, those two effects can clearly 
be distinguished. For very fast particles, this is not possible any more. This is particularly bothersome given that a 
pulse of gravitational waves will predominately excite w-modes. If we have a wave signal from a particle that grazes 
a neutron star with almost the speed of light, and we find that, indeed, there are some traces of a w-mode, it would 
be very difficult to state that the excitation of w-modes is a "real" signal and not an artifact due to the inappropriate 
initial data. 

To obtain an approximate answer, we again turn to the flat space case. In this limit, the equations governing the 
evolution of S and T reduce to two simple coupled wave equations with a source term for each that takes into account 
the presence of the particle: 

^ T + 4 ^ - 8. -I - RitW;^ . (35) 

Herein, we have omitted the terms proportional to L, since we let the particle move on a radial trajectory with 
constant velocity v, hence it is L = and 

R{t) = Ra+vt, (36) 

with Rq being the initial location of the particle. Furthermore, the normalized energy E of the particle is just given 
by the Lorentz factor 

It is interesting to note that the wave equation for S is totally decoupled from the one for T. However the solutions 
of (p3) and (pSl) have to satisfy the flat space Hamiltonian constraint which reads 



gj.2 J.2 dr \ 2 J r 



We now seek an exact solution of (34) that obeys the right boundary conditions at the origin and at infinity. Once 



found, we may use (38) to numerically compute the appropriate T. We state that a solution for (M) is given by the 



following series ansatz: 



/ '21+1+1 Ir, I f\l-2l-2 \ 

S{t, r)=A[Y. a. ^ ^ ^^,^^,^3 e(i?o + 1 - r) + ^ 6. ^^L±^ e(.-i?o-.^)j , (39) 



{Ro + vt)^'+^+^ 
where Q is the Heaviside function which satisfies 



0(a;) = ^' . (40) 

II, a; > 

Continuity at r Rq + vt requires that 

E«« = E^- (41) 

1=0 1=0 

The overall amplitude will be determined by A, hence we deliberately may set 

5]a, = ^6, = 1. (42) 

i=0 i=0 
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The amplitude A and the coefficients Ui and bi can be found by plugging (|3^) into (^4|). For A we find 



16n fiE^v'Y;^ ^^g^ 



and the coefficients and bi are determined by the following recursion relations 



2 (2i + ; + 3)(2i + / + 4) 

»-+^=°-" 2(» + l)(2z + 2^ + 3) (4^) 

^ 2(^-2i-2)(;-2z-3) 
^'+^ = ^'- 2(. + l)(2.-2/ + l) ■ 

It is interesting to note that, while the series in never terminates, the series in bi always terminates because one of 
the two factors in the numerator will become zero for some i. The series in converges if and only if |w| < 1. 

In Fig. H, we show the initial data obtained from (39) for a particle that is located at Rq — 500 km with different 



initial velocities. Here, instead of S and T, we plot the quantities rS and T/r, since from ( pq ) it is clear that only 
those expressions allow direct comparisons. For w = it is rS* = but the amplitude of rS grows rapidly when the 
particle's velocity approaches the speed of light, whereas the peak of T/r slightly decreases. In the ultra-relativistic 
limit rS totally dominates over T/r. 

The solutions above for S and T are not valid any more if we consider the particle in curved spacetime because 
they would violate the constraints. However, we can still use (^) as a prescription for S and then use the curved 
space Hamiltonian constraint ( p6| ) to solve for T. Furthermore, we can compute K from K = dS/dt and then use the 
momentum constraint ( |29| ) to compute K^. As long as the particle does not have any angular momentum and is far 
away from the neutron star, the thus obtained initial values should be a good approximation for a boosted particle 
on a Schwarzschild background. But if the particle has a large angular momentum, there will be additional source 
terms in the evolution equations and our approximation should break down. However, we are mainly interested in 
trajectories which come very close to the neutron star, and, hence, have only small angular momentum. In this case 
the above prescription still yields good initial data. 



V. NUMERICAL RESULTS 



To compare the code with known results, we first consider a particle in a circular orbit around the neutron star 
and compute the radiated energy at infinity. Numerically, this can be accomplished by evaluating the Zerilli function 
at some large distance. In terms of the metric variables S and T, the Zerilli function reads 

'-"^ ^2rT'+'Ji^lSl±I(^T-2r^s) . (46) 



l{l + 1) l{l + l)-2+^ \ ' ' r-2M 

The radiated power for some particular values of I and m can then be computed from 

dEim _ 1 {1 + 2)1 Zim. 2 
dt ~ 647r (/ - 2)! '"dT' 



(47) 



Cutler et al. have numerically computed the normalized gravitational power {M / fx)^ dEim / dt radiated by a 
particle orbiting a Schwarzschild black hole for various values of ?, m and Rq/M. In table II, they show the multipole 
components for Rq/M = 10. To compare the output of our code to those results, we chose the mass of the neutron 
star to be M = 1.99km and the orbit of the particle the be at Rq = 19.9km in order to obtain the ratio oi Rq/M = 10. 
The mass of the particle was set to ^ — 1 km. In Fig. ^ we show dZim/dt as a function of time extracted at r = 500 km 
for I = m = 2. After some wave bursts that come from the inappropriate choice of initial data, we see that the signal 
is periodic with a frequency of twice the orbital frequency of the particle. The amplitude is about 0.0076 which 
corresponds to a radiated power of (M//i)^ dE22/dt = 5.46 x 10~^ which is in excellent agreement with Cutler et al. 
who obtain a value of (M//x)^ dE22/dt = 5.388 x 10~^. The slightly higher value of our result may be due to the fact 
that at r = 500 km the Zerilli function is still somewhat off its asymptotic value which will be a little bit smaller. 
And, of course, a neutron star is not a black hole, since its oscillation spectrum is quite different. But as long as 
the particle does not excite the eigenmodes in a significant way, i.e. as long as the particle's orbital frequency is not 
in resonance with one of the eigenmodes, the emitted radiation should be quite similar to the black hole case. And, 
indeed, we find that we agree with all the even and odd parity modes compiled in table II of Cutler et al. within a 
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few percent, only for the case I = 3 and m = 1 we find the radiated power to be 5.9 x 10 instead of their value of 
5.71 X 10"^ 

If the orbital frequency is in resonance with one of the eigenmodes, the emitted radiation can drastically increase 
Ip^ . It is clear, however, that only the low-frequency modes such as the /-mode can be excited by this mechanism, 
the frequencies of the w-modes are much too high to be considerably excited. Besides, due to their strong damping, 
there is no sharp resonance for the w-modes at all, it is only for the weakly damped fluid modes that resonant mode 
excitation can occur. 

Therefore, the only way to possibly excite w-modes are very eccentric orbits, where the periastron is very close to 
the surface of the star or scattering processes with very small impact parameters. The investigation of the latter is 
the main objective of this work. Before presenting the results, we should point out some problems that are related 
with the use of an evolution code. 

Ideally, one should extract the waveform at r = oo which is clearly not possible on a finite grid. We have to record 
the waveform at some finite distance from the star. But this means that if the particle moves on an unbounded orbit 
it will eventually cross the location of the observer who will be located on a finite radius Vobs- Of course, the particle 
is always slower than the propagation speed of the gravitational waves, hence the observer will see the wave signal 
first before the particle crosses robs- However, when the particle is very fast and the observer is not far enough away 
from the neutron star, then the signal and the particle will cross the observer almost at the same time. That means 
that the signal will be a superposition of the "real" gravitational wave signal and the infiuence of the gravitational 
field of the particle itself. The farther the observer moves outwards the better the separation of the two components 
in the signal can be made. Furthermore, the amplitude of the signal remains constant for increasing r, whereas the 
influence of the particle decreases at least as 1/r. The effectiveness of this undesired influence of the particle strongly 
depends on the actual excitation strength of the neutron star oscillations. Smaller turning points of the particle 
will induce oscillations with higher amplitudes and thus make the influence of the particle undetectable. For large rt 
the induced fluid oscillations are so weak that they will totally be buried within the gravitational field of the particle. 
Figure || drastically shows the influence of the presence of the particle on the observed waveform. Here, we compare 
the signals that result for two scattering orbits with rt = 10 km and rt = 50 km. One can also see that the effect 
is much more pronounced in the waveform of T. But regardless of this unwanted feature it is still possible to infer 
qualitative statements from the obtained wave signals. 

Unfortunately, the influence of the source terms is even much worse for the Zerilli function, which should be used 
to present the results, since it is a gauge invariant quantity and directly related to the radiated power. This is not 
the case for our variables S and T. But since we are mainly interested in assessing which modes are excited by the 
scattering process, we do not have to compute the Zerilli function. It is enough to just examine the waveforms of S 
and T. 

In Figs. ^ through [ill, we show the waveforms of S and T for different initial radial velocities vq and different values 
of the radius of the turning point rt. In Figs, [l^ through |l^ we show the associated power spectra. 

For all cases considered, the waveforms always have the same overall features. The first part of the signal consists 
of an increase in the amplitude that results of the particle entering the strong gravitational field of the neutron star. 
After a few irregular oscillations which correspond to the particle curving around the neutron star, the final signal 
will be dominated by the quasi-normal ringing of the neutron star, which is confirmed by Fourier analyzing the signal. 

For slow particles the power spectrum shows one broad peak around the frequency that corresponds the the angular 
velocity of the particle at the turning point rt and some little sharp peaks that correspond to the fiuid modes of the 
neutron star (/-mode and lowest p-modes). From the waveforms it is clear that the excitations strengths of the modes 
decrease rapidly as the turning point rt of the particle moves farther away from the neutron star. 

When the particle reaches the speed of light, the picture changes drastically. For initial velocities vq < —0.6 the 
appearance of the w-mode becomes evident. In the ultra-relativistic case of vq — —0.97 (Fig. pT| ) which corresponds 
to E Ri 4.5 it is clear that the lion's share of the energy gets radiated through the first w-mode channel. Note, 
that in obtaining the power spectrum of Fig. |lj, we have cut off the initial radiation and have taken the Fourier 
transformation right where the w-mode ringing starts. It is interesting to see that even if the particle's turning point 
is quite far away (r^ — 50 km) the w-mode is still visible. 

All the calculations have been performed using a stellar model with central energy density of eo = 5.7 x 10^^ gcm~'^ 
which leads to a mass of 1.35 Mq and radius of 7.52 km. This model is right below the maximum mass that is allowed 
for the polytropic equation of state with F = 2 and k = 100 km'^. It is clear that the results will not change too much 
if we change the equation of state, or if we consider models with different masses. Only for very light neutron stars, 
the w-modes can become invisible again. For ultra-compact stars, we expect the long-lived trapped modes to be a 
dominant part of the power spectrum. 
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VI. CONCLUSIONS 



We have performed simulations of particles on circular orbits as well as on scattering trajectories. The circular 
orbits were done for code testing reasons, and we found that the radiated energies agree perfectly with those obtained 
by Cutler et al. who looked at a particle orbiting a black hole. 

For the scattering orbit, we discussed the numerical difficulties associated with an evolution code. However, the 
results are accurate enough to address the issue of the excitation of w- modes by orbiting particles. First of all, it 
became clear that in order to excite a significant amount of w-modes the particle's velocity at infinity must be 3> 0.5c. 
Secondly, the excitation strength of the modes rapidly decreases as the particle's turning point rt increases. This is 
in agreement with the results of Ferrari et al. [lO| , who report that no significant w-mode excitation is observable. 
However, they did not consider the case where the particle has an energy i? 3> 1, which is necessary if it were to 
excite w-modes. This was done by Andrade and Price for the odd-parity case ||l2[. Here, too, they find excitations of 
w-modes. 

The following question now arises: Can we infer anything about what happens in an astrophysical event? Our 
results show that the particle's initial velocity has to be incredibly high in order to produce a significant amount of 
w-mode excitation. It is clear that there do not exist any astrophysical events at all that could accelerate the particle 
(which, of course, represents some heavy extended object, like a planet or even another neutron star) up to some 
significant fraction of the speed of light. Hence, we might safely conclude that any astrophysical scenario that might 
be simulated by a particle orbiting a neutron star will not produce any detectable amount of w-modes. 

The alternative is then to excite w- modes from the collision of the particle with the neutron star. Our method 
could not be applied to this situation, but we could speculate about the possible implications. It has been suggested 
that in this case there might be some significant excitation of w-modes. In particular, simulations of a particle that 
spirals onto a ultra-compact constant density star where done by Borelli |^ who found that, indeed, both the trapped 
modes and the w-modes where excited, the latter much less, however, than the former. One serious drawback of 
those calculations, however, was that as soon as the particle had hit the surface of the star, the calculation was 
stopped. This lead to an overestimation of the high frequency components in the power spectrum which, of course, 
are represented by the w-modes. Since in this study, only odd parity modes were investigated, it is clear that the 
w-modes will always show up in the spectrum because they will not be screened by the presence of fluid modes what 
can happen in the even parity case. Now, even in the particle limit, it is not possible to simulate a collision since it 
not clear at all what happens if the particle hits the surface. 

A first step, though, might be to let the particle go right through the neutron star without being affected by the 
presence of the neutron star matter. If this creates strong w-mode excitations even for particles with low initial 
velocities, then one could conclude that, in a realistic scenario where the physics of the impact is included, there still 
might be w-modes present. Of course, in this case, it seems reasonable to infer that at least the fluid modes would 
be much more strongly excited than in a scattering event. 

Furthermore, our results seem to indicate that the presence of the w-modes in the signal is somehow related with 
the value of S. We have seen that for ultra-relativistic particles the values of rS greatly dominates the one of T/r. 
And it is only there where we flnd a signiflcant amount of w-modes. This particular role of S seems to be corroborated 
by the results of paper I, where it was shown that the occurrence of w-modes strongly depends on the chosen initial 
data. It was found that for conformally flat initial metric perturbations (i.e. = 0) the evolution shows practically 
no w-modes at all. However, if one takes the same initial fluid perturbation, sets Tq = and solves for S, suddenly 
one can observe a quite large excitation of the w-modes in the signal. This also agrees with results in [|j and [Q 

Still, the question remains, as to what happens in a realistic scenario. We have excluded the possibility of w- 
modes excitation by means of a realistic scattering process. However, it is not possible to relate our results to the 
merger process of a binary neutron star system since one cannot adequately described the system within perturbation 
theory. Only, if the final object does not immediately collapse to a black hole, it might wildly oscillate and emit some 
significant radiation through u;-modes. It will only be through nonlinear evolutions that this issue might finally be 
settled. 
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APPENDIX A: THE SOURCE TERMS OF THE PARTICLE 



The particle's energy-momentum tensor (|]) will be expanded using the orthonormal set of tensor harmonics 
{{yLU}^^^ ^^ given by ZeriUi 



oo ; 10 



~ X/ X/ X/ ^'^A[yWn]fii' 



1=0 m=-l A=l 

The coefficients can be obtained by means of the orthogonality relation 
where the asterisk denotes complex conjugation and 

{yLr = v''^v'"' [yLu- 



(Al) 



(A2) 



(A3) 



Because of the use of the inverse Minkowsky metric rjf^'^ to raise the indices, the inner product ( [A^ ) is not positive 
definite and it is sa = —1 for the harmoni cs w ith nonzero Oj-components {A = 2, 3, 4). 
By using the orthonormality condition (A2) we can compute the coefficients Pj^' through 



iA= I [ytr^T.^'^in 



(A4) 



We thus obtain the following set ( |Ts| with some corrections): 



[tt] 
[tr] 

' te 

tcj) 



llm 2z^ 

La — c 



2i^ 



4A 



re 

r(j) 



„2A 



r^2l{l 
2^ 


+ 1) 




r^2l{l 


+ 1) 


fi dr f' 




) 5{r^ 


'^'dt \ 


dr J 


2^ 




drdR 



r^2l{l 
2ifi 


+ 1) 


dt dr 
drdR 



I de d o^Aly* 
sine dr a$ ^'"^ dr dO 



Im 



llm 

to — 



ry^2l{l + l) dt dr 

fi dr 
^2l{l + - + 2)'dt 



5{r ~ R{t))—Y;^ 
5{r~R{t)) sme — — 



Tim 



llm 
'■10 



/i dr 
V^'dt 



Sir -Kit)) 



de 
d7 



Sir-Rit)) 



sin^ e 



dr de sin e dr 

2 



Y, 



de 



d^ 

d7 



dr 

2 



— 1 - sin^ e I — 



w,* 



y: 



dr 



^2l{l + l){l-l)il + 2) dt 



S{r - R{t)) sine 



d^ 



1 



sm 



de 



^e\dT 



with 



Xlr. 
Wlr^ 



2 1 -cote I —Yi„ 



(A5) 
(A6) 



11 



Here, Y";^, W*^ and X*^ are functions of the particle's angular position Q and $ parametrized by the coordinate time 
t, therefore all derivatives with respect to proper time r are to be understood as 

d dt d n\ ^ d , . „. 

— = — — = e^^E— . A7) 
dr drdt dt ^ ' 

Since the evolution equations are expanded using the Regge- Wheeler harmonics, which do not form an orthonormal 
set, we have to convert the above coefficients into the coefficients tj^* that would result from an expansion of the 
energy-momentum tensor into Regge- Wheeler harmonics. Using the relationship between the different sets of tensor 
harmonics given in JTsj l and | [T9|] we find: 

Ara jthn 

Zl — li 

Am ^ -f^/m 

2 ~ /7j''2 



V2 



Am ^ llr\ 

— , ^=l^A 



j/m _ Tin 

H — ^5 

An 
^6 



Am ' lln 



Am Am 

7 — — / 7 



vni + 1) 

,2 

^2Z(/ + l)(?-l)(/ + 2)''* 



2r.2 

A m Thn 



Am _ I flm, I / ^(^ + 1) 

^ "Til' +Va-i)a + 2) ' 



^/m ^Im 

y/2l{l + - 1){1 + 2) '° 



We now restrict the motion of the particle to the equatorial plane = f . In this case it is 4^ = and sin 8 = 1, 



and we can use the geodesic equations ( |3a[ ) and (|3c| ) to substitute all expressions containing derivatives with respect 
to proper time r. Furthermore, we can obtain quite simple relations for the derivatives of YJ*„: 



d_ 



YCr, = -miY:^ (A8) 
Y;^^{m'-l{l+l))Y:^. (A9) 
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This gives us a somewhat simpler set of coefficients: 

tt = e"^^5{r-R{t))Y:^ 
t\^ = -e^^^Vr5{r^R{t))Y{;, 



_ ^,LHl{l + l)-2m^) _ 

2im^iL^ d 



= ~^^\^Eiii + m-i){i + 2)^^'' " ^(*))ae^'" 

where = ^ is the radial velocity of the particle. The field equations also require the computation of the trace 

l,m 

with 

i'" = -e^^ t'™ + if - ^-^^^ 4" + ^ 4" ■ (All) 
Using the explicit forms of the coefficients we obtain 

t'"" = ^^i^" RW) e''-E + e''^^ YC^ , (A12) 

and by making use of the geodesic equation ( pb| ) we can reduce this expression to 

t'"' = ~e'-'-^6{r-Rit))Y;,^. (A13) 
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FIG. 1. Profile of the initial values for both kind of initial data. The left panel shows when To = 0. The right panel shows 
To when So = 0. Note, that S exhibits a discontinuity at the particle's location whereas T is continuous. 
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FIG. 2. Evolution of S. The left panel shows the evolution for To = 0, in the right panel So ~ 0. In the first case shows a 
burst of gravitational waves that propagate both in and outwards. The ingoing pulse gets reflected at the origin and travels 
back outwards again. In the other case, a wave originating from the interior of the neutron star travels outwards. Note the 
differences in the amplitudes of both cases. 

T T 




FIG. .'5. Evolution of T. The left panel shows the evolution for To = 0, in the right panel So = 0. In the first case, wc can 
see the built up of T which is disturbed by the reflected part of the wave that was sent out by the particle. The second case 
clearly shows that T basically has its "right" shape right from the beginning. 




FIG. 4. Evaluation of the Hamiltonian constraint during the evolution. Again, the left panel shows the evolution for To = 0, 
in the right panel So = 0. The Gaussian shape of the particle is clearly visible and has been chosen to be very broad for a 
better visualization. The particle is initially at rest and starts falling towards the neutron star. 
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FIG. 5. Plot of the variables S and T at the end of the evolution of the different initial data. The artificial radiation of 
the initial data has been radiated away and both variables have assumed their "right" values that is independent of the initial 
data. The difference in S is due to the fact that the initial data with To = contain much more radiation which excites the 
neutron stars to pulsations which in turn disturb the profile of S. The variable T is mostly unaffected by this. 
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FIG. 6. Initial data for a particle with different initial radial velocities v where 5* is given by ( p9[ ) and T is obtained by 
solving (pq). For small velocities, T/r dominates rS whereas for ultra-relativistic velocities, rS dominates T/r. 
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FIG. 7. Evolution of the real and imaginary parts of Z22 at 
r = SOOkm. After the initial radiation bursts the waveforms show 
periodical oscillations with twice the orbital frequency. The am- 
plitude is about 0.0076 which corresponds to a radiation power of 

{M/ij,f dE22/dt = 5.46 x 10"^ 
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FIG. 8. shows the difference of excitation strengths for the two different turning radii rt = 10km (solid line) and n = 50km 
(dashed line). The waveform of S (left panel) for rt = 10km is totally unaffected by the presence of the particle whereas for 
rt — 50km, the gravitational field shifts the amplitude of the signal to higher values. For T this effect is much more pronounced 
and can already be detected for rt = 10km. For rt = 50km, the oscillations of the neutron star are totally buried in the 
gravitational field of the particle. At t = 10ms the particle crosses the observer who is located at robs = 1000km. 
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FIG. 9. Waveforms of 5 and T for the three different turning radii rt = 10km, rt = 20km and rt = 30km. The initial velocity 
of the particle is vo = — 0.1 
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FIG. 10. Waveforms of S and T for the three different turning radii rt = 10km, rt = 20km and rt = 30km. The initial 
velocity of the particle is vo = —0.5 
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FIG. 11. Waveforms of S and T for the three different turning radii rt = 10km, rt = 20km and rt = 50km. The initial 
velocity of the particle is vo = —0.97 
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FIG. 12. Fourier transform of the waveforms for vo = —0.1 
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FIG. 13. Fourier transform of the waveforms for vo = —0.5 
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